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Abstract 

Modifications of the 7r° 27 decay amphtude by thermal effects 
have already been considered by several authors, leading to quite dif- 
ferent results. I consider in this paper the triangle diagram connecting 
a neutral pion to two photons in a constituent quark model, within 
the real-time formulation of thermal field theory and study the zero 
external momentum limit of this diagram. It appears that this limit 
is not unique and depends strongly on the kinematical configuration 
of the external particles. This non-uniqueness is shown to explain the 
contradiction between existing results. I end with some considerations 
suggesting that this decay amplitude may be significantly modified by 
the resummation of hard thermal loops, due to infrared singularities. 
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1 Introduction 



During the past two years, a lot of work has been devoted to the study of the 
relationship between the axial anomaly and the 7r° 27 decay rate at finite 
temperature, most notably by Pisarski, Tytgat and Trueman |||, |2|, |3|, |^, ^, ^ . 
The purpose of this series of papers was to explain the following basic fact: 
the coefficient of the axial anomaly is independent of the temperature while 
the amplitude for the it" — > 27 decay is modified. The problem was therefore 
to explain why the relationship that relates at zero temperature the pion 
decay amplitude to the axial anomaly ceased to be valid in a hot medium. 

This work has been initiated by a calculation of the pion decay rate in 
a constituent quark model, performed by Pisarski in the imaginary time 
formalism |jl], ^ . More precisely, it consists in the calculation of the triangle 
diagram connecting the pseudo-scalar to the two emitted photons, via a 
quark loop. This diagram is considered in the limit of vanishing external 
momenta. The result found in is that this diagram is proportional to 
m/T"^ where m is the mass of the quark in the loop and T the temperature 
of the heat bath, while the result found at zero temperature for the same 
diagram is proportional to 1/m. The consequence of this result is that the 
pion decay rate into two photons vanishes if the chiral symmetry is restored 
at high temperature, since m — > 0. 

The same diagram has been calculated in the real time formalism by 
[0, and also by Gupta and Nayak (GN in the following) in P who studied 
the zero momentum limit of this diagram. GN's result for this diagram in 
the zero external momentum limit is proportional to m/mT. The dramatic 
difference is the behavior of this decay amplitude as a function of the quark 
mass, because this behavior was crucial in Pisarski's calculation |jl], || to 
derive his conclusion about the pion decay rate in a hot chirally symmetric 
phase. 

The purpose of the present paper is to reconsider the calculation of the 
triangle diagram already studied by Pisarski and Gupta &: Nayak, in order 
to explain the discrepancy between the results they found. To that effect, we 
perform this calculation in the "retarded-advanced" version of the real time 
formalism, but we stay at a more general level than 0, 0, H i concerning the 
kinematical configuration of the external particles. In particular, we don't 
assume that external particles are on- shell. Like 01, we are interested in 
the zero momentum limit for this diagram. We arrive at the conclusion that 
this discrepancy is due to the non-uniqueness of the zero momentum limit of 
the considered Green's function. It appears indeed that this limit depends on 
the kinematical configuration of the external legs and that Pisarski and GN's 
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calculations correspond to very different configurations, GN's configuration 
being the most appropriate for the decay of a pion into real photons. Then, 
we come back to Pisarski's result about the pion decay rate in a hot chiral 
phase and show that, because of infrared singularities, it may remain valid 
in GN's kinematical configuration despite a different dependence in the mass 
m if one considers the correction provided by hard thermal loops. 

In section ^ we derive the expression for the triangle diagram in the 
retarded-advanced formalism, and its relationship with the pion decay rate. 
Then, we prove the existence of a limit of zero external momentum, in a 
sense to be made precise later. 

In section ^, we first give an expression for the zero momentum limit 
showing clearly that this limit is not unique and depends on the kinematical 
configuration of the external particles. The remaining of this section is 
devoted to the detailed study of this limit in three particular configurations. 
The first configuration studied corresponds to a situation where both of the 
emitted photons have zero energy: the zero momentum limit reproduces in 
this case Pisarski's result. The second important case is obtained with real 
photons and a pion at rest in the frame of the plasma: this case reproduces 
GN's result. Finally, a third simple case corresponds to the decay of a static 
pion into two static photons. 

In section ^, we study the implications of the above results for Pisarski's 
assertion concerning the pion decay amplitude in a hot chiral phase. Despite 
the fact that this assertion seems incorrect at first sight if one considers the 
physical situation in which the photons are real, the interplay of infrared 
singularities in this calculation makes the resummation of hard thermal loops 
necessary. The consequence of this resummation is to change the parameter 
playing the role of an infrared regulator. This has the effect of making the 
pion decay amplitude vanish in a hot chirally symmetric phase, even when 
one is considering the decay into real photons. 

Technical details are relegated to three appendices. In appendix^, we re- 
mind the reader of the potentially dangerous effect of changing the variables 
in divergent expressions since this is of some relevance for our calculation. 
Appendix ^ gives the general expression of the functions A{Ki,K2) and 
B(Ki, K2) that appear in intermediate steps of the calculations. Finally, 
appendix |^ gives some details about a few integrals that appear in this 
paper. 

2 Triangle diagram in the "R/A" formalism 
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2.1 71° decay rate 

The decay rate of pions in a thermal bath is related to the 7r°7r° retarded 
self-energy via the relation 



dN 



which gives the number of 7r° decays per unit time and per unit volume of 
the plasma, in the four momentum range dqod^q. The imaginary part of 
the 7r°7r° two-point function is a sum over all the possible cuts through the 
corresponding diagram, which means that this formula gives the total decay 
rate, i.e. the sum of the contribution of all the channels. In order to select 
a particular channel, one must look at the appropriate cut. 

Like H, ^, I use a linear sigma model (see |10| for instance) where the 
fermion fields are constituent quarks, in which the mesons are coupled to 
quark fields as indicated by the following Lagrangian 



(2) 



I consider two flavors of quarks and = 3 colors. The t matrices are 
normalized with to = 1/2 and Tr{tath) = Sab/"^- This coupling is invariant 
under the chiral symmetry SU{2)^ x SU{2)^ . When this symmetry is spon- 
taneously broken, the a field acquires a non vanishing vacuum expectation 
value|^ (a), which gives a mass m = g {a) to the constituent quarks. In this 
model, the decay of pions in two photons appear only in the discontinuity of 
the three loop 7r°7r° self-energy. Indeed, each external pseudo-scalars must 




Figure 1: Self-energy of the pseudo-scalar involved in the decay in 27. 

be connected to a quark loop, and these two loops must be linked by the two 
photons. Then, among all the possible cuts, one must consider the cut that 

^This vacuum expectation value can be identified with the pion decay constant for 
two flavors at zero temperature. At nonzero temperature, they differ somehow. Anyway, 
both of them vanish when the chiral symmetry is restored. 
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crosses the photon propagators (see figure |lj). Making use of the cutting 
rules for the "R/A" formalism |jll|, we find that the cut depicted on figure 
|l| contributes: 

_1 f (fKi f d^K2 , 
~2J J2^J (2^ 

X {27r f 6{Q + K,+ K^) (Q, K,, K2)t''^^ ^'^ {Q, K.^K^) , (3) 



'.2^e{k'i)5{Kl)2T,e{kl)5{Kl 



^ARR , 



where T [Q, Ki, K2) is the triangle diagram connecting the pseudo-scalar 
to two photons. This object will be the subject of our study from now 

RAA 

on (r^^, is closely related to the previous one). In fact, two diagrams 
contribute to this one-loop 3-point function because of the possibility of 
crossing the photons in the final state, as outlined on the figure 0. In order 




+ 




Figure 2: 1-loop triangle diagrams contributing to tt° 
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to take the two configurations into account, it is sufficient to calculate in 
detail the first one, and then add the term obtained by interchanging the 
indices (1,^) ^ (2,i^). 



2.2 Matrix element 



Let us first give the value of the vertex function T . Using the Feynman's 
rules established for the "R/A" formalism (see p^ , [l^), a straightforward 
calculation gives for one flavor of electric charge e: 



AmN e^ge^yap kfk^ 
X / ■0^{nAp" + k"2)s\P-Ki)s\p) Disc S''{P + K2) 

{p°)S^ {P - K{)S'' {P + K2) Disc S"" (P) 
+n^{p" - kDS'' {P)S'' {P + K2) Disc5''(P -Ki)] 
+ {Ki,fi)^{K2,iy) , (4) 
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where m is the mass of the quark running in the loopQand S (P) = i/{P^- 
m? lb ip°0"^) the scalar part of the retarded (advanced) quark propagator. 
In the fohowing, we can forget about the retarded or advanced labels for the 
denominators. Indeed, to recover the correct prescriptions, it is sufficient to 
perform at the very end of the calculation the substitutions: 



^3 



k°n - 2«0^ 



(5) 



It is worth recalling that the discontinuity of the quark propagator gen- 
erates a Dirac's delta function 



Disc S (P) = 2tt e{p°) 5{P^ - m^) 



(6) 



which enables us to do easily one of the integrations. To be more definite, it 
is convenient to use these Dirac's functions to perform the integration over 
the variable so that we are left with a three-dimensional integration: 



r;ir(K3,i^i,i^2) 



4mA^ (?g e^^yaji kfk^ 



p+k2 



2P-K2 + Ki 2P ■ {Ki + K2) + Ki - Kl 
i i 



-2P ■ Ki + Kl 2P-K2 + Kl 



+ 



-e{-e) 



2wp„fe, 2P-Ki- Kl 2P ■{Ki+K2)+ - Kf 



(7) 



where we denote ujp = ^/(p^ + m^). This expression^ of the vertex function 
will be the basis of further considerations. 



2.3 Existence of a zero external momenta limit 

We are interested now in the zero momentum limit of this vertex function 
in order to understand the origin of the discrepancy between Pisarski's and 

^Because the vertex coupling the pion to the quark loop is g^^, the result is proportional 
to the mass of the quark. If we were in a chirally symmetric model (m = 0), the Dirac's 
trace would be vanishing. 

^The reader who may wonder why we don't replace P + K2 by P in the first term and 
P — Kl by P in the third one is referred to appendix M. 
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GN's result. Only two of the three external momenta are independent ones 
due to the energy-momentum conservation: therefore we choose to consider 
Ki and K2 as independent momenta and replace everywhere^ by —Ki — 
K2. A priori, taking the limit Ki,K2 — > is a very intricate task since we 
have to take to zero the eight components of these four-vectors. In order 
to simplify without reducing significantly the generality of the result^, we 
will assume that the size of the eight components is controlled by some scale 
A, and this parameter will be the only one taken to zero. This amounts to 
write: 

= XKi , K2 = XK2 , (8) 

where the components of Ki^2 are fixed and of order unity. By this sub- 
stitution, we are lead to considering the limit when A — > of a univariate 
function F{\), the Ki^2 playing the role of constant parameters. 

We now want to show that the integral appearing in Eq. (0) has a finite 
limit when A — > 0. If we recall Eq. (|^), we can see that this integral is the 
sum of six terms (three terms, plus the terms obtained in the symmetrization 
with respect to the external photons), each term behaving like in the 
limit A ^ 0. Therefore, in order to obtain a finite result, we must expand the 
integrand in Eq. (^) up to the order A'', and show that we have cancellations 
among the various terms in order to eliminate the orders A~^ and A~^. 

The order A~^ is easy to obtain, since we can drop the A dependence in 
the statistical functions to extract it, which gives: 



2P ■K2 2P- {Ki + K2) 



+ 



2P ■Ki2P- {Ki + K2) 



2P ■Ki2P- K2 

+ ^ {k2,u) = Q. (9) 



As we can see, the cancellation of the order A ^ is in fact a consequence of 
the energy-momentum conservation (it works because we have replaced K-^ 
by -Ki - K2). 

^From now on, we drop the explicit reference to the argument in F^^ . 

^The most general case would be a situation where the eight components of 7^1,2 are 
arbitrary functions of the parameter A, vanishing when A ^ 0. For our purpose, it is 
sufficient to restrict ourselves to a linear dependence in A of these components. 
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The cancellation of the order A~ is a consequence of the parity properties 
in Ki and K2 of the vertex function. Indeed, looking at Eq. (0), it is rather 
straightforward to check the identity: 

r;r(i^i,K2) = r;r(-^i,-i^2) . (10) 

Making use of the variable A, it can be rewritten as: 

r;r(A,i^i,i^2) = r;r(-A,i^i,i^2) . (11) 

In other words, the one-loop vertex function is an even function of A. This 
implies automatically that the terms of order A~^ in the Laurent's expansion 
of the integral are vanishing. For this cancellation to occur, it is essential to 
perform the symmetrization with respect to the external photons. 
Therefore, if we write: 

(A,i^i,K2) = AmNe^ge^^^pk'^k^ f^f (A, Ki, i^s) , (12) 
then limA_>o r^^^('^; -^1; ^2) is finite. 



3 Non-uniqueness of the limit 

3.1 Generalities 

After some tedious expansions^, we find: 



lim 



111/ 



(2^ 



3 1 - 2np{u)p) 



B{k^,k2) <(cup) 



1 



TT ^ ^ 



(13) 



where we denote V± = {ujp, ±p). The functions A and B are quite intricate; 
since their detailed expression is not really helpful here, they have been 



®At this stage, once we have proven the existence of the hmit A — > 0, we can speed up 
the calculations by making use of some computer algebra system Uke Maple for instance. 
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quoted in the appendix The fact that the above expression still depends 
on Ki and K2 means that the value of the zero momentum limit depends 
upon the path chosen to reach the point Ki = K2 = in momentum space. 

The non-uniqueness of the zero momentum limit in this case should not 
be a surprise. Examples of such a phenomenon are weh known in thermal 

RA 

field theory. For instance, the same calculation applied to the IIqq compo- 
nent of the photon polarization tensor in massless QED leads to 

hm Ilt^{\,K) = 4e2 / n'ip) " . , (14) 

which is nothing but the HTL contribution to this function. Here also, the 
residual dependence upon K indicates the non-uniqueness of the limit. In 
both cases, this remaining dependence on how the small momentum limit 
is reached implies that the corresponding term in an effective Lagrangian is 
non-local. 

There is though an important difference between Eq. ( [T3| ) and the HTL 
amplitudes. The hard thermal loop approximation consists in retaining only 
two orders in the expansion in powers of A (the lowest order is trivially van- 
ishing due to momentum conservation). In the case of Eq. (|l3|), we have 
combined two diagrams so that the second order is also vanishing. We 
therefore need to calculate the third order of this expansion, and this is why 
the functions A and B are much more involved than what is usually encoun- 
tered in hard thermal loops (Eq. (14) for instance). As a consequence, one 



may expect that the effective 7r°77 coupling near the critical point exhibits 
a nonlocality of a completely different nature]]. 

Before going on with some specific kinematical configurations, a com- 
ment is worth concerning the zero temperature limit of Eq. (|l3|). Since for 
m > we have lim^^g (wp) = lim^^g n'^ (wp) = limy^g n'^ (wp) = 0, the 
zero temperature limit is trivial: 

Jim_^^r;r(A,i^i,i:2) = ^mN^ge^,o.pK4 j • (^5) 

As one can see, the integral is now totally independent of the kinematical 
configuration of the external particles. Therefore, the fact that the numerical 



^Let us recall that near T = Q, the nonlocality of the anomalous couplings is found to be 
closely related to that of hard thermal loops ^, ^ . More precisely, HTL- like amplitudes 
are encountered in thermal corrections at the order / fi in a low temperature expansion. 
Near the chiral phase transition, we are in the opposite limit T 3> /vr, and it is likely that 
new nonlocal terms appear. 
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coefficient in front of the zero momentum limit of this diagram may not be 
uniquely defined is a purely thermal effect. 



3.2 Space-like photons 

A first possibility is to consider the situation where /c°2 = while ki^2 7^ o. 
This corresponds to external space-like photons. In this particular case, the 
functions A and B become much simpler: 

A{Ki,K2) = -3(p • ki)\p ■ k2)\p ■ ik + k2)f 
B{K,,k2) = -iP ■ kifip ■ k^fip ■ ik + k)? ■ (16) 

Plugging these expressions into Eq. ([T^), we find: 

lim ^"IT (A, i^i , i^2) = 4mNe^g e^^af} Kk'^ J -^^^ 
/3 l-2n^(a;p) , 3 n'^jujp) 1 n'^juj,) \ 

We can perform at this point the analytic continuation of Eq. (|5|). Since 
the functions A{Ki,K2) and B{Ki, K2) exactly cancel the denominators of 



Eq. (13), this analytic continuation does not introduce any imaginary part 
in the result. This fact is a consequence of a result proven by Evans ||T5|] , 
according to which all the retarded/advanced Green's functions are equal if 
the external energies are set to zero. 

The angular integration is trivial here since it just amounts to multiply 
the result by Air. It remains to perform the integral over p = \\p\\- This 
integral cannot be performed analytically if m 7^ 0, but we can consider 
performing an expansion of the result in powers of m/T, assuming m <^T. 
In fact, replacing m by zero in the expression inside the brackets, we can see 
that the integral over p is infrared-safe without the need of this mass. As a 
consequence, the first term of the expansion in powers of m/T is trivial to 
extract: 

+00 

dp 



hm T^^^''{X,Ki,K2)='imNe^ge^,^pk'^k^ f 
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Integrating by parts in order to get rid of the inverse powers of p, we obtain: 

lim (A, KuK2) = - ^^^^^ e^yap fc^ 

+00 

j dxHx)h";{x)(l + o(^y^ , (19) 


where we denote np{x) = l/(exp(x) + 1) and x = p/T. Making use of 

n'^'(x) = 6nl (x) - 12n^ (x) + Tn^ (x) - (x) , (20) 
and of Eq. (|40| ) in appendix we finally findQ: 

hmr^, (A, Ki,K2)= ^1^2 (^1 + (^^ j j , (21) 

which is equivalent to formula (11) of Q. Therefore, we have shown that 
Pisarski's result, obtained in the imaginary time formalism with external 
momenta set to zero right from the beginning, corresponds in fact to a zero- 
momentum limit taken with space-like external photons. 

This fact can be interpreted as follows: since in the imaginary time 
formalism the energy component of four vectors is a discrete quantity, the 
only possible way of taking the "zero momentum limit" in this formalism 
is to first set the external "energies" to the discrete value zero, and then 
consider the limit of zero three momenta. The above analysis shows that 
the limit is unique once the external energies are set to zero (the dependence 
on ki^2 lias disappeared in Eq. (p!7|)), which implies that the imaginary time 
formalism leads to a uniquely defined limit that coincides with the result 
obtained here with space-like photons. 

A remark is worth concerning the paper [^] by Baier, Dirks and Kober, 
who reproduced the result of |Q] in a somewhat different framework. In- 
stead of calculating the triangle diagram in a particular model, they consid- 



ered the Wess-Zumino-Witten [18, 19] functional near the chiral symmetry 



restoration. Intermediate steps of their work involve the calculation in the 
imaginary time formalism of a function where the external momenta are 
set to zero. It seems that this technical analogy with [|l| is the reason 
of the agreement. Since the zero momentum limit of the 7r°77 triangle is 



*The formula ( |4C| ) of appendix ^ naturally leads to the quantities C(~2) and (^'( — 2). 
In order to simplify the result, we use the identities C(~2) — and (^'(—2) — — ^(3)/47r^ 
(see for instance jLql). 
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not uniquely defined, a complete calculation of the Wess-Zumino-Witten 
Lagrangian near the chirally symmetric phase should be extremely careful 
when using the imaginary time formalism (or avoid it), in order to get the 
correct nonlocality for the couplings contained in this functional. 

3.3 Real photons 

Gupta and Nayak choosed to consider the decay of a massive pion at rest 
in the frame of the plasma into two real photons. This choice corresponds 
to the constraints ^1 + ^2 = and kf = k2 = — ||fci,2||) implying some 
simplifications for the functions A and B: 

A{k^,k2) = 16fc?V[(a;,fc?)2 - {p ■ k^)^f 

Biki,k2) = ^ktWiP ■ kif[iu;,klf - ip ■ fci)2] , (22) 



and for the vertex function: 



hm r;^^ (A, ki,k2)= 4mNe^g e^,^p k'^k^ f 

A — ^0 J 



d?p 



(27r)3 

3 1-2n^(cup) 1 <(u;p) 1 ^(a;^) {p ■ 



,(23) 



where we denote k = 1 1 the unit vector in the direction of the emission 

of the first photon. The analytic continuation of Eq. (^) generates a term 
d{ujp — {p ■ k)'^). Anyway, since ujp > p, the Green's function T^^^ remains 
real|^. 

As one can see now, the angular integral is not defined if the quark mass 
is vanishing, due to a collinear singularity. This could have been expected 
since we are looking at the emission of real photons. The angular integration 
gives the expression: 




2ujj 2a;3 \^ 2p \ujp - p 



+ ^T-f^ .(24) 



This resuh is very similar to the situation encountered in the calculation of hard 
thermal loops which have an imaginary part only if some external momentum is space- 
like. Here, it still works for external lines on the light-cone thanks to the mass m in the 
loop. 
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Besides the potential collinear singularity, another dramatic difference of 
this case with respect to the previous one lies in the infrared behavior of 
the integral. It is now impossible to take the limit m ^ in the expression 
inside the brackets because the integral over p would behave like dp/p^ at 
small p. This means that the expansion of the integral in powers of m/T 
begins with a term in 1 /mT, to be compared with the 1 /T^ of the previous 
situation. Integrating by parts the above equation, we can transform it into: 

hmr,. (A, K,,K,) = ^vaNe ge.^^pk, k^.j ^^^^^ ^^^^ In (^-— j , 

(25) 

which is equivalent to the result given by Gupta and Nayak for the decay 
of a static pion into two real photons (see formula (2.12) of The first 
term of the expansion in powers oi m/T is: 

limrr(A,K„A-,) = 9^.,„„,kffc?(l + 0(^)) ^ (26) 
3.4 Photons at rest 

Another simple case is the situation where the emitted photons are both 
massive and produced at rest in the frame of the plasma (they can subse- 
quently decay into lepton pairs). Therefore, the kinematical constraints we 
must enforce are now fci,2 = o while 7^ 0. With these constraints, the 
functions A and B become trivial: 

A(i^i,it2) = 0, i?(i^i,i^2) = 0, (27) 

so that we have: 

limr;r(A,i^i,J^2) = AmNe'ge,^^,k?k^,J ^ - ^M^p) . (28) 

Again, the angular integration is trivial and for the remaining integral on the 
variable p we can only perform an expansion in powers of m/T. The analytic 
continuation of Eq. (^) has no effect on this result. Here also, this integral is 
infrared divergent if we put m = in the integrand. As a consequence, the 
result of the integration behaves as 1/mT instead of 1/T^. More precisely, 
we have: 

limrr(A./?.K.) = ^.„„„,*f*,f(l + 0(^l)) , (29) 
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3.5 Epilogue 

The above particular examples have demonstrated clearly the non-unique- 
ness of the zero momentum limit of the triangle diagram responsible for the 
pion decay in two photons. Moreover, the particularity of the first situation 
must be emphasized: when one expands the integral in powers of m/T, there 
is a cancellation of the terms of order 1 / mT so that the first non vanishing 
terms is of order A closer look at the functions A and B in appendix 

^ and at Eq. ([l^) indicates that the point where both = and = 
is quite exceptional, because the functions A and B vanish faster at small p 
when /c° = = 0, implying that the integral is not infrared sensitive. When 
at least one photon has k° ^ 0, then at least two powers of p are replaced 
by m in the small p behavior of A and 5, so that the expansion starts at 
the order 1/mT. Therefore, generically, the 7r°77 effective coupling does not 
vanish in the limit of chiral symmetry restoration m ^ 0. 

4 IR sensitivity of tt^ 27 and hard thermal loops 

4.1 Preliminaries 

The behavior of the decay rate of the 7r° into 27 when the chiral symmetry 
is restored is closely related to the behavior of the F^^, function in the 
limit where the mass m goes to zero. The above study shows how this 
behavior depends on the kinematical configuration of the external photons. 
In particular, we observe that the imaginary time calculation performed with 
external momenta set to zero does not correspond to the physical situation 
where the emitted photons are real, but rather to a situation where the 
photons are both space-like. The fact that the imaginary time calculation 
does not correspond to real photons could have been expected thanks to the 
absence of any collinear singularity in this approach. 

The problem is now that GN's situation, which seems more physical be- 
cause the photons are assumed to be real, leads to a very different behavior 
for the triangle diagram at small m. Indeed, Pisarski's result behaves like 
m/T^ and therefore vanish in the limit of chiral symmetry restoration. On 
the contrary, GN's result behaves like ra/mT and therefore tends to a non 
vanishing constant when we consider the same limit. The question is there- 
fore: is the conclusion that the 7r° — > 27 decay rate vanishes if the chiral 
symmetry is restored at finite temperature correct, since it has been derived 
using the result for space- like photons ? At first sight, it seems that this 
conclusion is erroneous, because it makes more sense to consider the result 
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established for real photons in this context. 

4.2 Infrared sensitivity and hard thermal loops 

Nevertheless, another aspect of the problem is to be considered, which may 
have important consequences in the limit of chiral symmetry restoration. 
Indeed, as seen above, the zero momentum limit in the case of real photons 
contains a strong infrared divergence, which gives the factor 1/mT (instead 
of once regularized by the mass m. This means that the integral over 

the loop momentum is dominated by the soft scale. More accurately, the 
momentum p, of order m, becomes softer and softer as one approaches closer 
to the chiral symmetry restoration. Therefore, there is a point when the 
loop momentum is soft enough to justify the resummation of hard thermal 
loops 1^, ^] on the quark propagators, as outlined on figure ^. 

Since we have two coupling constants e and g in our model, we can define 




Figure 3: 7r°77 Green's function at one loop in the HTL expansion. 

two soft scales eT and gT. But since the coupling constant g is related to 
strong interaction while e comes from electro-magnetic interactions of the 
quarks, one may expect that loop corrections involving the constant g are 
the dominant ones. As a consequence, we will consider only loop corrections 




Figure 4: Dominant topologies contributing to the HTLs of figure ^. 
involving the cr or tt, as shown in figure ^. Looking at the Lagrangian in 
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Eq. (^), we see that the couphng of the a field to the quark field is —ig 
while the coupling of the tt to the quarks is (77^ (I don't write here the 
isospin matrices since they appear in such a way that the end result for the 
quark self-energy is proportional to the identity in flavor space). Moreover, 
in order to derive the contribution of a and tt loops to the HTL correction to 
the quark propagator, we can neglect the mass m of the constituent quarks 
since m <C gT when we approach the critical temperature. We obtain for 
the retarded self-energy at HTL order: 

it (P) - . 2 / MO + MOl 2 

-^T-RA^^JiTV = / TT^Tv^ ^, „ f 7 PI , (30) 



1^ ^7 (27r)3 21 P-L 

where we denote L = (1,Z). As one can see, the sum of the above two con- 
tributions is precisely equal to the standard result of QED with replaced 
by g^. As a consequence, we know already all the properties of the effective 
propagator obtained by the resummation of the above self-energies: this re- 
summation introduces a cut-off of order gT in loop involving these effective 
propagators. 

For the sake of completeness, we can give expressions for the HTL part 
of the above vertices. Starting with the 'jqq vertex, we find 

r755- (o P-P-O)'- -iea^ f Kil)+nAl)] W .3^. 

' i (2vr)3 41P.LR.L ^ ^ 

for the contribution of the a field, and 

' J {2^r AlP.LR.L ^ 

for the contribution of the tt field, where we denote R = P + Q. If we add 
the two contributions, we can see that it is equal to the QED HTL vertex 
with two factors of e replaced by g. Exactly in the same way, we can obtain 
the HTL contribution to the vertex TT°qq. In fact, the result is obtained by 
substituting in Eqs. (^) and (^) the last power of e by (7 and the matrix 
7'^ by ^7^. Since 7^ anti-commutes with the Dirac's matrices, we see that 
the product of matrices entering in this vertex is always proportional to 
fjj^ = = 0. Therefore, the 7r°gg vertex does not have a HTL contribution 
at the scale gT (it may have one at the much smaller scale eT though). 
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4.3 Effect of HTLs on vr" 27 

This resummation has the effect of giving a thermal mass to the quark, 
and this thermal mass remains constant]^ as we approach the point of chiral 
symmetry restoration (i.e. limy_^j,^ rrij, ~ gTc), while the constituent quark 
mass m goes to zero (limj,^^^^ m = 0). In other words, the thermal mass 
will become the relevant infrared regulator as soon as m <C rrij. {i.e. the 
relevant infrared regulator is always the biggest one available). Moreover, 
the thermal mass has the property of not modifying the chiral properties 
of the propagator which means that if m 0, the thermal mass will not 
change anything to the fact that the Dirac's trace vanishes. On the basis 
of these arguments, one can expect a modification of the m/mT behavior 
in the case of real photons. Indeed, if we track the origin of the various m 
factors in this result, we see that the mass m in the numerator comes from 
the Dirac's trace, and is closely related to the fact that the chiral symmetry 
is broken by the mass m. On the contrary, the thermal mass rrij, does not 
break chiral symmetry. Therefore, this factor m at the numerator remains 
unmodified by the resummation of the thermal mass. The mass factor in the 
denominator has a completely different origin: it comes from the infrared 
sensitivity of the integration over the loop momentum. This infrared scale is 
affected by the resummation of the thermal mass. As a consequence, we may 
expect that the mass m is replaced by the thermal mass in the denominator 
but not in the numerator, when m <C tti^. 

Therefore the resummation of the thermal mass m^, would lead to the 
behavior mlm^T for real photons, near the chiral symmetry restoration. 
As a consequence, the result according to which the pion decay rate for the 
channel 7r° 27 vanishes in the chiral phase at finite temperature survives. 
More precisely, the situation that emerges from the constituent quark model 
used in this paper is summarized on the figure ^. Assuming that the chiral 
symmetry restoration is a second order phase transition, the constituent 
quark mass is a function m = m(T) of temperature that vanishes at a certain 
critical temperature Tc- We can also define two additional temperature 
scales T\ and T2 for which m(Ti) = T\ and m{T2) = 5T2, respectively. 

Then, for T G [0,ri], the temperature has negligeable effects, and the 
decay amplitude is very close to the zero temperature one: it behaves as 1 /m. 
When the temperature reaches values of order of Ti, thermal corrections 

^° Assuming a second order phase transition, we expect the coupling constant g to depend 
only logarithmically upon temperature, while the vacuum expectation value of the sigma 
field, responsible for the mass m of the constituent quarks, vanishes as a power of T — Tc 
at the critical temperature Tc- 
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Figure 5: Shape of the temperature dependence of the n° 2j decay amplitude in a 
constituent quark model. The dotted curve (1) is an extrapolation of the T = result with 
a T dependent quark mass. The dotted curve (2) is the result obtained with real external 
photons when one takes into account thermal corrections, but not the resummation of 
hard thermal loops. 

become important, which has the effect of replacing the factor l/m? by 
1/mT. As a consequence, the decay ampHtude behaves Hkc 1/T for T € 
[Ti,T2]. Finally, in the domain [r2,Tc], the resummation of thermal masses 
plays the dominant role in the regulation of infrared singularities, and the 
decay amplitude eventually vanishes for T = Tc since it goes like m/gT^. 

5 Conclusions 

First, we have seen that the difference between the results of Pisarski and GN 
can be interpreted as an effect of the non-uniqueness of the zero-momentum 
limit of the triangle diagram at finite temperature. It appeared also that 
Pisarski's result, originally derived in the imaginary-time formalism, corre- 
sponds in fact to a situation where the emitted photons are space-like. The 
most physical situation corresponding to the case where the emitted photons 
are both real, a superficial analysis tends to invalidate the result according 
to which the 7r° 27 decay rate vanishes in a hot chirally symmetric phase. 
Nevertheless, this calculation seems incomplete in the chiral limit since 
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the loop integral is now sensitive to soft momenta, which means that the 
resummation of hard thermal loops may have important effects. Taking 
them into account will change the infrared regulator, and modify GN's result 
in such a way that it now vanishes in the hot chiral phase. 

I will add also a word of caution concerning the imaginary time for- 
malism. Since the zero momentum limit of thermal Green's functions is 
not unique and depend on the location of the external legs with respect 
to the light-cone, this limit cannot be handled correctly in the imaginary 
time formalism (there is no light-cone in an Euclidean formalism). Indeed, 
in such situations, this formalism gives a number which corresponds to one 
particular way of taking the limit, but which is not necessarily the most 
appropriate for the problem under study. Moreover, the information about 
the nonlocality of the corresponding effective coupling is lost. 

Among the related topics that seem worth studying, I would mention the 
case of the box diagram appearing in 7r°a — > 77, since one can expect here 
also to have a nonlocal effective coupling in a hot chirally symmetric phase. 
Another interesting aspect is to find a work-around for the limitations of 
the imaginary time formalism in the derivation of the anomalous processes 
near the phase transition from the Wess-Zumino-Witten Lagrangian. 
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A Side-effects of hazardous translations 

A superficial inspection of the expression of the vertex function given in 
Eq. (0) indicates that it may be advantageous to change the integration 
variable in the first and third term. Namely, changing P + K2 into P in the 
first term and P — Ki into P in the third term would allow the factorization 
of a common statistical weight (n^(u;p) — 9{—e))/2ujp. Besides this factor- 
ization, the major advantage of such a transformation would be to eliminate 
the K12 dependence (i.e. the A dependence) in this statistical factor: as a 
consequence, the subsequent expansions in powers of A would become much 
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simpler. This technique can be apphed in GN's situation where it leads 
directly to Eq. (|2^), without the need of performing a cumbersome integra- 
tion by parts (the reason for that lies in the fact that no derivatives of 
are generated in this second approach, because the statistical prefactor does 
not contain the expansion parameter A). This method works also for the 
third case studied in section ^ Nevertheless, we avoided its use to derive 
the general limit presented in Eq. (|l3|) because this transformation is not 
always legitimate. 

The reason why the above transformations are sometimes illegitimate 
comes in fact from the infrared sector. To be more definite, let me focus on 
the situation of Pisarski, since for this one the above changes of variables 
are not allowed. The particularity of this configuration is that k° = k2 = 0. 
Therefore, the denominators in Eq. (0) are combinations of 2p - ki — k\ and 
2p • ^2 + ^2- As a consequence, the expansion of these factors in powers of A 
generates powers of kf /{p ■ ki). A simple counting shows that the order A'' 
behaves like / dp/p^ in the infrared region, even if we keep m strictly positive. 
Therefore, each individual term in Eq. (Q) is strongly infrared divergent at 
p = Q. Of course, a conspiracy of the three terms occur in order to cancel 
this divergence so that the final result is finite. If one performs different 
translations on the three terms, then these finite terms are modified. The 
correct answer can only be obtained when the same transformation (or no 
transformation at all) is applied to the three terms. 

An alternative way to see that this transformation is not legitimate in 
Pisarski's situation is as follows. We have seen on the GN's case that an 
integration by parts relates the expression obtained without this transfor- 
mation to the expression one would obtain by making use of it. Therefore, 
let us integrate by parts the result of Eq. ([TtI) : 

2n^(wp) ^ ^n'^{ujp) 




uj^ 2 



-oo 








2^ 3 In^ iMin) 6_M_4 

+ 00 

2n,(m)^ I ^^n,(^,)--n,(m)_ ^^^^ 



^m? J 2p^ujjj 



As one can see, some terms in np{m) appear in this integration by parts, 
which are absolutely necessary to ensure the infrared finiteness of the in- 
tegral. There is no way of transforming this integral in order to have the 
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temperature dependence only in a factor 1 — 2np{ujp), as it would be if the 
changes of variables P + K2 P,P — Ki — > P were possible]^. 



B Functions A{Ki, K2) and B{Ki, K2] 



A{Ki,K2 



l.o3 T,o5 (n 1.0 
1 2 \ 1 



^2) (^1 ^2 "I" ^2 ^1) ' 



fcf (p • kif + 2 fcf kf (p • fci) (p • fcs) 
fc2 (p . kif + 10 fcf kf (p • fci) (p • k2) 



-2kr kr {p-kif-kfkf (p-fca) 



.o3 T o2 
„o T o2 7,2 , 



+2 kl kf ki {p-ki){p- fca) -^kl kf (p • fci)^ kl 
+6 fcj' fcf fcf (p • + 4 A;f A;^ kl (p • fci) (p • fcz) 
-10 yfcf kl {p-kif + 8 kf k^ (p • ki) (p • fca) 
-2 A:f A:^ fef (p • fci) (p • fca) + 10 kf k^ kf (p • kif 



-3 A:f k^ (p • fea)' - 3 kf k^ kf (p • fca) 
-lOAf fc2 (p.fc2)' + 5A;f fc2 (p-fca)' 



to, 
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4A:?^ {p-k^)ki{p-k2t 

-Skfkl {p-ki){p-k2f -{p-ki 
-2 fcf kl (p • fci)^ (p • fca) - 20 Af I 
-3 fe^^ fc? (p • kif - 26 /tf kf (p • fc; 
-9 k°2 kf (p • fca)^ + 7 Af Af (p • kif (p • fca 



2(p.fci 
3(p.fci 



(p-fc2 



4 ;.o3 



-36 kf kf (p • kiY (p • fca) + 3 A^ (p • k2Y kf (p • k 
-12 kfkf {p-kif 
+kfkf {p-k,f{p-k2f 



18 kf 



(p-kifk 

<o2 7 o3 



4 kfkf (p-fci)(p-fc2) 



-12 kf ip ■ kif kl (p • k2f + 3 kf kl (p • k2f 



-4 fe° kl (p • fci)'' (p • fca) + 2 A? (p • kif ki{p-k 
+7 kf kl (p • kif kl-bk°2 kf kl (p • kif 



^'^The condition to be able to eliminate the np{m) terms is 
ujp d p 



1™ a 

p->0 p op 



dnpB{Ki,K2)Y[- 



which is satisfied in GN's case but not in Pisarski's one. 



0, 



(34) 
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-3 kf kl kl (p ■ kif (p ■ kif - 5 kf kl kl {p ■ kif 

-kfkl{p-kif + 2kf{p-kifkl 

+7kl (p • kxf kf kl-6k^ kf kj (p ■ fci) (p • fea)^ 

-9 k^ kf (p ■ kif kl (p • k2f -4k^kl kf (p ■ kif (p • fes) 

-4 fc^ fcf (p • kif kl (p • k2) - 12 A;2° fcf kj (p • fci)' (p • k2f 

+(p ■ fci) kf [6 (p • kif kj kl (p ■ k2f 



+3 (P • fcl) fe^ kl (p • fc2)" - 4 ^2° (P • fcl] 



klkl{p-k2) 



+ Q{p-kifkt + ^{p-ki)''{p-k2)kl 



2 7.o3 



+2(p-fei)^fc?(p-fe2rA;^ 



+ 17 (p ■ kl) kf kf (p • k2T + 20 (p • fci)* A;r (p • k2) kl 
+ 4 (p ■ fei)' kf kl (p • fc2)' + 56 (p • ^2^2 ^o2 . ^^y 
+ 13 (p • fci)^ fc^' kf + 13 fef (p • kif kf 
+ 34 kf (p • fci)^ (p • k2) kf + 2 kl kf (p • ^2)^ 
+ 31 k"2 (p • kl) kf (p • fcs)' + 40 k"2 (p • ki)^ 
+ 16 (p • fci 



(P-A;2)' 



kfkl {p-k2f + Qkfkl {p-k2f 



+ 42 kf (p • fci)^ (p • k2f kf + 14 (p • fci) kf ki (p • A;2) 



^2 7.o2 



o2 



+ 11 (p • fci) kf (p • fc2)' fci - kf (p ■ fci)-^ kl 

- 2 ^-2°' (p • fc? (p ■ k2f + 4 ^2°' (P • kif kl (p • fe2)' 

+ 4 A;2"2 (p • ki)^ kl (p • k2) - k°2 (p • kif kl kl 

+ 4 kl (p • fci)' fci (p • fc2) fe? - 2 ^2° (P • fcl)' kl (p • fc2) fc? 

- 4 fcf (p • kl)" kl (p • fes) + 2 A;2° (p • A:? fcf (p • fea)' 
+ 8A;§(p-fei)2A;?fc?(p-fc2)^^ 



~{p-ki)'kl Wkf{p.ki)'{p.k2) 



22kf{p-ki)\p-k2 



+ 24 kf (p • (p • fc2)' + 3 kf (p • fci)^ + 8 kf (p • (p . A;2)' 



+ 23 kf (p • fci)" (p • k2)' kl + 44 kf (p • fci)^ fc^ (p • fc2; 
+ 2 kf (p • fci) fc? (p • fc2)' + 19 kf (p ■ fci)' A;? (p • fc2)' 
+ 12 kl (p • fci (p • ^2)=^ - 4 kl (p • fci)^ fc? (p • fc2)' 
+ 25 kl kf (p • A;2)^ + 76 fc^ (p • kl) kf (p • fc2)^ 
+ 2 fc^ (p • ki)^ kl (p • fc2) - 3 kl (p • fci)' fc? (p • fc2)^ 
+ 8 fc^ (p • fei)^ (p • fcs)' fei - Af (p • kl)" (p • fc2)' fc? 
+48 kl (p • fci)' (p • fc2)^ fcf + 6 A;^ (p • ki)^ kl (p • fes)^ 
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+4 ip ■ fci) k1 kj {p ■ k2f + K{p- k^f kl 
+{p ■ kif kl kl [p ■ k2f - 2 (p • fci) kl {p ■ k2f kl toi 

+6{{p ■ k2) + {p ■ fci)) k"2 iP ■ fci)' iP ■ k2f [k"2 {P ■ kif 

+3 (p • kif k"2 {p ■ k2) 

+ 4 k"2 iP ■ kl) {p ■ k2f + 2 ^2 (P ■ k2f + 2kl{p- k2f] u^l 

- 3 (p • kit {p ■ k2f {{p ■ k2) + {p ■ fcl)) (3 {p ■ fci) + ip ■ k2)) 
+ {kl,ki)^{k"2,k2) , 



B{Ki,K2) = kl kf [kl ip ■ k2f + k°{p- kif 

-kl {p ■ kif [{kl + k"2) ip ■ k2f - kl {p ■ ki){p ■ (fci + k2)) 

-{p.kif{p-k2f{p-{ki + k2)) 
+ {kl,ki)^{kl,k2) . 



C Calculation of /o+°° x"" ln(x)(exp(x) + 1)-Pdx 



In the section we need to evaluate integrals of the form: 

'n,p = 



ln(x) 

In. . = I dx ^ ^ , (35) 







{e^ + 1)^ ' 



where n,p are positive integers. The starting point is to expand (e^ + 1) ^ 
in powers of e~^, which gives: 



/I 

n,p 



p-1 

(v- 1)! ^ 

V-f^ -^^ m=l 1=0 



( — 1)P +°° 

I ) ^(-l)™^ap_i,im* / dxx"ln(x)e- 



(36) 



where the numbers are the coefficients of the polynomial 

p-1 

Qp-i{x) = {x - l){x - 2) ■ ■ ■ {x - p + 1) = ttp-i^i X* 

i=0 



We need then 



+ 00 



dx x" ln(x) e 



An - n!(7 + ln(m)) 



m 



n+l 



(37) 



(38) 
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where 7 is the Euler's constant, and An are integers defined recursively^ by 

Ao = An+i = n! + (n + l)An . (39) 

It is now straightforward to collect the various pieces in order to obtain the 
following expression 



(2*~" - l)C'(n + 1 - i) - 2^-" ln(2)C(n + 1 - i) 



+(2--_i) (^-7) Cin + l-i) 



(40) 
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